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Abstract 


A new formulation of helicopter rotor thick- 
ness, noise for hover and forward flight, is dis- 
cussed. The parameters required for this formula- 
tion are rotor motion, planform and airfoil thick- 
ness distribution. A computer program has been 
developed to calculate the pressure signature due 
to blade thickness for a helicopter in arbitrary 
motion. Comparison with high-speed helicopter 
tests shows good agreement with calculations when 
the observer is in or near the horizontal plane 
in which the rotor disc lies. Characteristics 
of thickness noise are illustrated by numerical 
examples indicating strongly that the high-speed 
blade slap may be due primarily to the thickness 
effect. The methods of Deming and Arnoldi are 
discussed as the special cases of this technique. 


I. Introduction 


An unexplained phenomenon of rotor noise 
generation is the high-speed blade slap. This 
noise which appears at high advancing tip speeds 
has several distinctive characteristics. It is 
impulsive, that is, its signature has a high crest 
factor, and is directional towards the forward 
region of the flight. The peak pressure is in 
the rotor plane. The sound intensity appears to 
increase substantially as the relative tip speed 
approaches sonic speed in the medium. This noise 
is sometimes referred to as the compressibility 
noise for the obvious reason that the compressi- 
bility effects become apparent at high tip speeds. 

It is the purpose of this paper to suggest the 
thickness noise as the possible mechanism of high- 
speed blade slap. The theory developed by 
Farassat IT) has been used to study the character- 
istics of the thickness noise. Some comparison 
with experimental measurements is made and reported 
here. The agreement with experiment and 
the characteristic features of the thickness noise 
suggest this mechanism as the source of the blade 
slap. 


II. Thickness Noise Theory 

The present theory is based on the solution 
of the wave equation: 
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where p is the acoustic pressure, x is the 
source time, and v^ is the normal velocity 
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of the blade surface which is defined by the 
equation f(y, x) =0. The density of the undis- 
turbed medium is p Q and the speed of sound is 
denoted by c. In Equation (l), 6(f) is the 

Dirac delta function. This equation has several 
forms of solution which were originally studied by 
Ffowcs Williams and Hawkings One form based 

on a collapsing sphere which was also proposed by 
Ffowcs Williams and Hawkings and developed by 
Farassat (3), seems to be appropriate for numerical 
calculations. The derivation of this solution is 
given in (l). The solution of Equation (l) used 
here is 


. / p cv 

i+7Tp(x,t) = -r— I — d T d.T (2) 

dt J r sin Q 

where x and t are the observer positions and 
time, respectively, v n is the normal velocity of 
the blade surface and r = 1 3c — y"l where jr is 
the source location on the blade. In Equation (2), 
F is the curve of intersection of the blade sur- 
face and the surface g = x- t+ ^-=0 which for 
a fixed x and t is a sphere whose radius 
collapses at speed of sound c. The angle 0 is 
between the outward normal to the body and the 
direction r = 3c - ‘y . Figure 1 illustrates the 
definition of some of the parameters. 



Figure 1. Illustration of the Blade System Inter- 
secting the Collapsing Sphere g = 0 
Forming T-curve . 


The normal velocity v of the airfoil can be 
written in the terms of helicopter speed and rotor 
angular velocity SI. Consider a rotating system 
of axes r) such that R^-axis is parallel to the 
chord of a blade, the positive direction being from 
the leading edge to the trailing edge. The rotor 
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system is in the n-j_n 2 ”P lane "with the origin at 
the center of rotation. Then 

F f (t"h ) 

Vn ”yi+F'2(n 1 ) (n 2 Q - V + V 3 n 3 (3) 

where FCn^) is the thickness function of the 
airfoil and V and are the helicopter speed 

along n 1 - and ru -axes > respectively. Here n 
is the unit outward normal to the blade surface. 
One can approximate the F-curve with the inter- 
section of the sphere g = 0 and the mean air- 
foil surface. The upper and lower surfaces of the 
airfoil contribute equally to the thickness noise 
so that one obtains, from Equations (2) and (3), 
the following: 


p(x,t) 
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The summation notation is used here to denote 
the integration over all the blades of the rotor 
system for which the intersection with g = 0 
occurs. Note that the axial velocity V 3 can- 
cels out in the integrand in Equation (U) when 
combining the integrals over the upper and lower 
surfaces of the blades. 


The condition 6=0 occurs when the surface 
of the blade and the sphere g = 0 are tangent. 

It may be shown that 00 Equation (2) can be 
written in the following form: 


l+irp(x,t) 
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where A = [l + Mn^ - 2M n cos 0] , M n = v n /c 

and Z is the surface that the T-curves form 
in the space for fixed x and t. The Z-surface 
is given analytically by f(y, t - £.) = 0. It is 
the locus of the points on the blades whose signal 
arrive simultaneously to the observer at the time 
t. From the relations 

A 2 = (1 - M n ) 2 + 2M n (1 - cos 6 ), 

1 + M 2 - 2M cos 0 £0, 
n n 

one can see that the only possible singularity 
of Equation (5) is when M = 1 and 0=0 simul- 
taneously. This condition requires more careful 
consideration ( ^ ). In case ^ 1, we see that 
the condition is a removable singularity of 
Equations (2) and (1+) so that a very small por- 
tion of Z-surface near such a point can be 
neglected. This can be incorporated in the 
numerical scheme and, thus. Equation (H) can be 
used for calculation of thickness noise. For 
the case M = 1 and 0=0 simultaneously at 
a point on the Z-surface, either smaller mesh 
must be taken during numerical calculation or the 
contribution from the region near this point must 
be found analytically and be included in the 
numerical scheme. The differentiation with res- 
pect to the observer time can be performed 
numerically. 


The extent of the Z-surface or more pre- 
cisely the extent of that region of the Z-surface 
with large source strength, the observer position 
and the observer time interval of interest as 
well as the source time scale determine the com- 
pactness of a sound generator. 

Apart from the neglect of the differences in 
the retarded time, the compactness character of 
the sound source makes it possible to relate the 
acoustic pressure to the global parameters of the 
source such as the net force on the body or the 
net rate of mass injection by the motion of the 
body surface. Thus, the surface structure of 
the source does not enter the acoustic calcula- 
tions directly, but the structure indirectly in- 
fluences these calculations by affecting the 
global parameters. For a compact source the 
Z-surface and the body surface should be nearly 
identical and the differences between the retarded 
time from the observer should be much less than 
the typical time scale of the source. For non- 
compact sources, the surface structure does enter 
the acoustic calculations (3), as is also evident 
from Equations (2), (h) , and ( 5 ). From this dis- 
cussion, one can conclude that even for moderate 
tip Mach numbers, with the observer in or close 
to the plane where the rotor disc lies, the com- 
pactness assumption is not justified and one 
should use Equation (4). 


III. Previous Work 


The earliest published work on the thickness 
noise is by Ernsthausen (5), He measured the near 
field sound radiation from a high-speed model 
propeller with symmetrical airfoil section and 
zero pitch angle. He explained the shape of the 
acoustic pressure signature by sweeping the 
pressure pattern around the airfoil past the ob- 
server. He discovered several features of the 
thickness noise such as the increase of the sound 
intensity and appearance of higher harmonics at 
higher tip speeds. 

Deming developed a theory for static pro- 
pellers ( 6 L. In this case, because of the peri- 
odicity of the signal, the spectrum of the sound 
can be obtained easily. His equation and the 
boundary condition is equivalent to the time 
Fourier transform of the single Equation (l). His 
result is, therefore, equivalent to the Fourier 
transform of Equation (U) for the case of a static 
propeller or a hovering rotor. Because of unavail- 
ability of high-speed computers , Deming was ob- 
liged to make further simplifying assumptions to 
get analytic expressions for his computations. 
Recently, Hawkings and Lowson (7) published a 
work which includes a thickness term identical 
to Deming 's using a slightly different approach. 

Using the work of Billing^ Arnoldi ob- 
tained a theory for moving and static propellers 9). 
His theory is too detailed to be discussed here, 
but the bulk of manipulations is due to the com- 
plex helical pattern of the source motion. The 
basic idea behind the theory was explained by 
Hawkings and Lowson (10) using the modern mathe- 
matical tool of generalized functions. Arnoldi' s 
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solution is identical to the solution of the 
equat ion 

JL JLE- _ = aL 

c 2 3 t 2 3T ^i 


-\Po U i V o 6( ^B (T)) } 


where is the source speed, which is at the pos- 

ition y^( T ) and has the volume V . In other 
words, Arnoldi's work is restricted to compact 
sources only. Also, his analysis is restricted 
to uniform axial motion of the propeller which 
cannot be applied to helicopter rotors. In time 
domain, one can easily solve the above equation 
numerically. 
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Figure 2. Theoretical variation of acoustic 
pressure signature due to thickness with blade tip 
speed for a hovering helicopter. Ho. of blades = ; 
Rotor diameter = 10 m. Hub diameter = 1.1* m. Blade 
aspect ratio = 10.75, Blade chord = 0.1* m, airfoil 
type-biconvex parabolic, observer in the plane of 
the rotor 50 meters from the center, of rotation, 

= tip Mach number (tip velocity /c). Thickness 
ratio = 10 percent . 


IV. numerical Calculations and Comparison 
With Experiment 


50 r 


Some of the characteristic features of the 
thickness noise are discussed here by numerical 
examples. Also, the results of two helicopter 
tests are compared with calculations. The compu- 
ter program with some discussion of the theory 
will be published at a later date. The basic 
features of the present method is that all the 
calculations are performed in the time domain by 
constructing the T-curves for each source time 
T and performing the integral in Equation (4) 
by finite difference method. 

For Figures 2 and 3, the rotor system is 10 
meters in diameter with two rectangular blades. 

The airfoil section is a biconvex parabolic arc. 
The computer program is not restricted to airfoil 
sections with sharp leading edge as Figures 5 and 
6 show. The chord length is 0.4 meter and 
the airfoil thickness ratio is 10 percent. The 
aspect ratio of each blade is 10.75. 

Figure 2 shows the effect of changing tip 
speed from tip Mach number Mt (tip velocity/c) 
of 0. It to 1.1 for a hovering helicopter. The 
acoustic pressure signatures at a point 50 meters 
from the center of the rotation and in the plane 
of the disc are given for one period of the 
signal. One should note the considerable in- 
crease in amplitude with increase of the tip speed. 
Starting with approximately equal positive and 
negative peaks, the negative peak becomes domi- 
nant from about M. = 0.6. For sonic and low 
supersonic tip Mach numbers, there are two large 
positive and one large negative peak. The nega- 
tive peak for in the range 0.4 to 1.1 follows 

closely a variation given by exp (8.9** M-fc 2 ). The 

appearance of higher harmonics is also evident in 
Figure 2. For higher tip Mach numbers, one ob- 
tains a signal which is concentrated in a shorter 
portion of the period with larger peaks. 

Figure 3 indicates the change in the pressure 
signature with observer on a sphere 50 meters in 
radius and changing elevation from the plane of 
the disc up to 45° elevation in 15° increments. 

The tip Mach number is 0-9 and the helicopter is 
in hovering position. The signature becomes less 
impulsive with the increase in elevation and the 
amplitude decreases considerably indicating 
strong directionality in the plane of the rotation. 

Figure 4 corresponds to the condition of 
Figure 2(f), but with rotor hub diameter of 8.4 
meters (84 percent of the rotor diameter) result- 
ing in blades with aspect ratio of 2 instead of 
10.75. This figure shows that it is the tip 
region of the blades which is responsible for 
most of the noise generation. 

Figure 5 is for a two-bladed rotor system of 
a helicopter at forward speed of 87.5 m/sec (170 
kts). The rotor diameter is 13.42 meters, the 
blade chord is 0.686 meter and the rotor RPM is 
324. The airfoil thickness ratio is 9-3 percent. 
Figures 5(a)? ( "b ) , and (c) show the theoretical 
pressure signatures for blades using a biconvex 
parabolic arc , an HACA four-digit and a super- 
critical airfoil, respectively. At emission time, 
the observer is 80 meters ahead of the helicopter 
which is flying at 15 m altitude. These 
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Figure 3. Theoretical variation of acoustic 
pressure signature due to thickness with observer 
location for a hovering helicopter. Rotor para- 
meters as in Figure 2, observer 50 meters from the 
center of rotation, = 0.9. 


signatures are calculated to show the effect of the 
change in thickness distribution of blade section. 
It is seen that the biconvex airfoil is the least 
noisy and the supercritical airfoil the most noisy 
of the three sections. It must be mentioned, 
however, that the use of supercritical airfoils 
may make it possible to reduce the thickness ratio 
near the tip and, thus, result in both improved 
performance and better noise characteristics. 


OJ 



<D 

U 

0 

CO 

CO 

(L> 

u 

Cu 


50 

0 

; J 



1 period 

-50 

- 

Period = 50.6 msec 

-100 

- 


-150 

- 



Figure The influence of "blade tip region on 
the generation of thickness noise. Rotor hub 
diameter = 8.4 m, = 0.9. Other parameters 
as in Figure 2. Compare this figure with 
Figure 2(f ) . 
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Figure 6 is for the rotor system of Figure 5, 
but with an HACA four-digit airfoil section. The 
observer position is as for Figure 5- Figures 6(a) 
and (b) are for forward speed of 72.0 m/sec (l40 
kts) and 87.5 m/sec (170 kts), respectively. The 
corresponding experimental acoustic pressure signa- 
tures for a test helicopter with the same rotor 
parameters and operating conditions, obtained at 
NASA Langley Research Center, are superimposed on 
the figure. The experimental signatures also in- 
clude the tail rotor noise. It is seen that the 
thickness noise theory can explain the observed 
experimental signature. The discrepancy may be 
due to the errors of numerical differentiation in 
using Equation (it). It may also be due to the 
neglect of the contribution from the pressure 
sources on the blades which is not considered in 
this paper (7). For the observer in or close to 
the plane where the rotor disc lies (near or far 
field), thickness noise appears to be of the cor- 
rect order of magnitude for experimentally observed 
pressure levels of high-speed blade slap. 



experimental 



Figure 6. Comparison for theoretical and experi- 
mental acoustic presure signatures for a test 
helicopter. Airfoil type - NACA four-digit. Other 
parameters as in Figure 5. 


V. Conclusions 


Figure 5- Effect of change of thickness distri- 
bution on the theoretical acoustic pressure sig- 
natures. No. of blades = 2, Rotor diameter 13. **2 m 
Blade chord = 0.686 m, RPM = 324, Helicopter 
altitude = 15 m, observer location = 80 m (ahead) 
of helicopter) at emission time. Thickness ratio = 
9.3 percent. 


This paper presents a new formulation of the 
thickness noise theory which is not restricted to 
either compactness or far-field assumptions . Due 
to the complex geometry of the T-curves or the 
E-surfaces for a rotor system, it is difficult to 
draw general qualitative conclusions based on 
Equation (4), for the arbitrary motion of the 
helicopter. Numerical examples which are worked 
out indicate that: 
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1. The blade noise due to thickness at high 
tip speeds is impulsive in character and can have 
large amplitudes of the order observed for high- 
speed blade slap. 

2. The thickness noise is directional with 
peak pressure in the rotation plane towards the 
forward region of flight. 

3. The tip region of the blade contributes 
mainly to the generation of the thickness noise. 
The signature is directly proportional to the 
thickness ratio 

k. The variation in the thickness distribution 
of the blade section can have substantial in- 
fluence on the noise radiation. 

To achieve noise reduction, one may reduce the 
blade thickness ratio and change the thickness 
distribution near the tip region of the blades. 

The effect of the planform changes has not been 
studied here, but this offers another means of 
controlling the radiated noise. 

Based on the fact that the present theory 
exhibits the observed features of high-speed 
blade slap of helicopters, thickness noise is 
strongly suspected to be the source of the blade 
slap at high relative tip speeds. Additional 
investigation is needed to verify this conclusion. 

The nonlinear propagation effects which may be 
important if the propagation distance is large, 
has been studied by Hawkings and Lowson (7). 

The change in the wave form when included in the 
analysis, gives better agreement with experimental 
results according to these authors. 

The first author acknowledges the support from 
NASA Grant NGR-09-010-085, entitled "Noise Reduc- 
tion." The computer program was developed with 
the help of Mr. G. H. Mall of Computer Sciences 
Corporation. 
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